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Abstract 

Mixtures of cold bosonic atoms in optical lattices undergo phase separations 

on different length scales with increasing inter-species repulsion. As a general 

rule, the stronger the intra-species interactions, the shorter is this length 

scale. The wealth of phenomena is documented by illustrative examples on 

both superfluids and Mott-insulators. 
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The first experimental realization of the superfluid to Mott-insulator phase transition of 
cold bosonic atoms in optical lattices (OLs) [1], following the theoretical suggestion in [2], 
has boosted the community to study the physics of bosons in OLs, see, e.g., the recent review 
[3] and references therein. Recently, mixtures of bosonic species in OLs have attracted ample 
attention, see, e.g., Refs. [4-12]. Obviously, the physics of bosonic mixtures trapped in OLs 
is more intricate than that of the single-species case [1-3,13,14]. For instance, due to the 
inter-species interaction additional quantum phases appear [5-7,9]. 

Phase separation is a basic phenomenon distinguishable particles can undergo with in- 
crease mutual repulsion. How does a mixture of bosonic atoms undergo phase separation in 
an OL, instantly or, perhaps, step by step? Is there a difference between phase separation 
of two superfluids or of two Mott-insulators? In answering these questions, we show below 
that a bosonic mixture in OLs undergo a wealth of phase separations on different length 
scales. As a general rule, the stronger the intra-species interactions, the shorter the length 
scale in which phase separation in OLs can take place. More details are given below. Phase 
separation, by definition, happens in real space, which is 'where' we choose to attack this 
problem. 

Our starting point is the many-body Hamiltonian describing interacting iV = Na + N B 
bosons, Na bosons of type A and Nb bosons of type B, in an optical lattice, 

Hab = H a + H b + W ab . (1) 

Here Ha is the usual single-species Hamiltonian (for the A-type bosons), containing one- 
body Jia and two-body interaction terms. Similarity, H B describes the S-type bosons. 
Wab contains the inter-species two-body interaction terms. In what follows we employ the 
common contact interaction between all bosons in the mixture, see, e.g., Ref. [15], where the 
corresponding interaction strengths are denoted by qa, Qb and Qab and are proportional to 
the intra- and inter-species s-wave scattering lengths. 

As mentioned above, we are to attack the properties of a bosonic mixture in the OL in 
real space. Recently, a multi-orbital best-mean-field ansatz for single-species bosonic systems 
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has been derived [16] and led us to predict a wealth of quantum phases and excitations of 
strongly-interacting bosons in OLs [14] and novel phenomena associated with fragmentation 
and fermionization of bosons in traps [17]. Anticipating that an approach similar in spirit 
would be valuable for mixtures of bosons in OLs, we proceed to construct one. 

Consider the Hamiltonian (1) in the absence of intra-species interaction, Wab = 0. The 
two species are then independent and can be treated separately. The total wavefunction of 
the system ^ ab = ^a^b is just a product of the wavefunctions ^ a and ^> b of the individual 
species. Take species A and write a multi-orbital ansatz for ^ a- This amounts to attaching 
an orbital to each of the Na bosonic atoms. Generally, we may take rt\ A-type bosons to 
reside in one orbital, </>i(x), n 2 A- type bosons to reside in a second orbital, 4> 2 (x), and so 
on, distributing the Na atoms among 1 < n or b < Na orthonormal orbitals. More formally, 
we use the ansatz ^ A = «s{</>i(xi) • • • 0i(x ni )</> 2 (x ni+1 ) • • • </> 2 (x ni+n2 ) • • • ^^(xjvj}, where 
S is the symmetrization operator and Xj is the coordinate of the i-th A-tjpe boson. Taking 
the expectation value of the Hamiltonian H A with respect to leads to the single-species 
multi-orbital energy functional [16]: 

^ = e{"*/ 4>*(x)h A (x)Mx)dx + 9A ni{ni ~ 1} / |<Mx)| 4 dxj 

^orb p 

+ £ 2g A n i n j / |0 i (x)| 2 |0 J -(x)| 2 dx. (2) 

i<j 

In order to find the ground state of A bosons in the subspace of all possible configurations 
^a, one has to minimize Ea with respect to its arguments. These are the number n or b 
of orbitals in which the A-type bosons reside, the occupations {n,i\ of these orbitals and 
the orbitals {&} themselves which have to be determined self- consistently. Performing this 
minimization results in a system of n orb coupled, non-linear equations which is solved self- 
consistently [14,16]. Next, the S-type bosons are treated in exactly the same manner, and 
one only needs to introduce a separate notation for the B-tjpe orbitals, their number and 
occupations, and the B bosons' energy-functional which we respectively denote by {tp p }, 
rriorb and {m p }, and Eb- We point out that due to the distinguishability of the A- and 
S-type bosons, no relations are assumed between the </>'s and -0's orbitals which, therefore, 
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are a priori allowed to overlap. 

With all ingredients at hand, let us switch on the inter-species interaction Wab- We now 
employ \& ab = ^a^b as the ansatz for the wavefunction of the mixture including inter- 
species interaction, analogously to utilizing ^/a as the ansatz for interacting single-species 
bosons [16]. Taking the expectation value of the full Hamiltonian H AB with respect to ^ ab 
readily leads to the multi-orbital Bose-Bose energy functional: 



Finding the minimum of Eab amounts to minimize it with respect to the variational pa- 
rameters of Ea and of E B which were already discussed above. But now, the inter-species 
interaction Wab couples all orbitals which, as we shall see below, leads to very many effects. 
Before proceeding, it is gratifying to recognize that the two-component Gross-Pitaevskii ap- 
proach, see, e.g., Ref. [15], is a specific case of Eab ' ab) when all A-type orbitals are alike 
and all £>-type orbitals are alike, namely n or b = m or b = 1. 

The multi-orbital Bose-Bose energy functional (3) can be employed to study demixing 
in optical lattices by following the ground state as the inter-species repulsion grows. Let us 
examine some of its properties. Consider a given bosonic mixture with fixed intra-species 
interaction strengths gA, 9b and inter-species repulsion g A B (we shall only consider repulsive 
mixtures here). We can minimize the energy functional (3) and find the ground state, i.e., 
the set of orbitals {4>i}, {ip p } and their respective occupations {ni}, {m p } that minimize 
Eab- For g A B *C 9a,9b the single-species terms Ea and E B dominate Eab because the 
inter-species interaction term (the third term on the r.h.s. of Eq. (3)) is much smaller. Con- 
sequently, both species can completely mix and spread all over the optical lattice. Moreover, 
the physics of each species is determined by its own parameters almost independently of the 
other species state. Increasing g AB influences the physics of both species. In order to reduce 
the inter-species interaction energy in (3), the orbitals {&} and {ip p } reduce their overlap 
in space. When g A B gA, 9b we expect this overlap to vanish, which means that the two 
species completely separate and occupy different regions in space. In between these two 




p 



(3) 
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extreme cases, as we shall see below, there can be many intriguing possibilities of phase 
separation on smaller length scales. In fact, even when the two species fully separate in 
space there are differences, e.g., between two superffuids and two Mott-insulators. 

In a general bosonic mixture, the A and B bosons have different masses, filling factors 
and intra-species interactions, and they can possibly feel different optical potentials. Here, 
to capture and explore the richness of possibilities, still maintaining a coherent line of expo- 
sition, we consider both A- and £>-type bosons to have the same mass, the same intra-species 
repulsion gA,9B and filling factors, denoted hereafter by /a,/b, and to feel the same op- 
tical lattice. Additionally, we examine mixtures in one-dimensional optical lattices. It is 
convenient to introduce the dimensionless coordinate x in which the one-body Hamiltonian 
(of each species) reads h = + \E R cos 2 (x), where E R stands for the depth of the 

optical lattice in recoil energy units. Throughout this work we take E R = 20. The coupling 
constants Qa-,Qb and gAB are now related to the scattering lengths and the confining param- 
eters of the transverse directions [18]. In one-dimension one often expresses the strength of 
(intra-particle) interaction in terms of the dimensionless parameter 7^4 which is the ratio of 
interaction and kinetic energies. In the above units one has 7^ = ttj^ and similarly for the 
.B-type atoms. 

Let us begin by studying the phase separation in the simplest case, i.e., of two (weakly- 
interacting) superfluids. Taking gA = 9b = 10~ 3 and gAB = 10~ 4 and, for instance, Ja,b = \ 
we minimize the energy functional (3) and find that it is minimized when all A bosons reside 
in one orbital and all B bosons reside in another orbital, i.e., n or b = m or b = 1. Therefore, for 
these parameters the best mean-field of the bosonic mixtures is the two-component Gross- 
Pitaevskii approach [15]. How do the corresponding orbitals look like? They are delocalized 
over the whole lattice, see Fig. 1A. Obviously, no demixing occurs. Next we consider a 
much stronger repulsion, gAB = 1, and reminimize Eab- We again find that Eq. (3) is 
minimized with one orbital per species; the orbitals, however, now separate in space, see 
Fig. IB. This is the simplest phase separation possible. It is interesting to examine the 
interface region between the now separated species, see Fig. IB. It shows smaller densities 
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which grow and saturate as one moves away from the center of the interface. Clearly, by 
lowering the species' densities in the interface region the energy is minimized. Why can the 
density of each species at all be lowered? Because after phase separation had occurred each 
species remains a superfluid and is, thus, compressible. Repeating the above calculations 
for the same interaction strengths and the filling factors of /a,b = 1 we obtain the same 
scenario described above and in Figs. 1A,B. 

The above phase separation of two superfluids is on the longest possible length scale, 
i.e., of the order of the lattice size. We notice that the length scale available for phase 
separation is simply determined by the sizes of the orbitals in which the different species 
reside. This brings about the following idea. If we could make the sizes of these orbitals 
smaller then we might be able to realize states where the bosonic species are phase-separated 
on a shorter length scale. How can we make the size of the orbitals smaller? By increasing 
the intra-species interactions each species would independently fragment, namely it would 
occupy more orbitals of smaller sizes, see Ref. [17]. Along these lines, we consider two 
strongly-interacting superfluids, specifically we take qa = 9b = 0.1 with Qab = 0.1, and 
Ja,b = \- Minimizing E A b we find that: (1) each boson resides in a different orbital, (2) 
each orbital is spread over two lattice sites, and (3) orbitals of different species are alternating 
in position, see Fig. 2A. In other words, phase separation on a much shorter length scale 
has been achieved. Taking now a much stronger inter-species repulsion, g A s = 10, the 
strongly-interacting superfluids fully demix, see Fig. 2B. Next, we take the same interaction 
parameters but smaller filling factors, Ja,b — \i to allow for more space for each boson. 
Indeed, the size of each 'domain' is now three lattice sites, see Figs. 2C,D. Additionally, we 
observe that the spatial densities are maximal in the center of each 'domain' and minimal 
in the interface of two 'domains', in reminiscence of the finding for two weakly- interacting 
superfluids shown in Fig. IB. Decreasing the densities to Ja,b = where p is an integer, 
it is anticipated that larger 'domains' whose sizes are p lattice sites can be formed. 

So far, we considered phase separations on length scales which are larger than the size 
of one unit cell. Let us move towards phase separations on smaller sizes. Consider the case 
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of /a,b = \ and take — Qb — 0.1. By minimizing Eab f° r 9ab = we obtained that 
all bosons reside in different orbitals and that each orbital is spread over two lattice sites. 
Of course, this is a system of two independent species, each being a strongly-interacting 
superfluid. What happens as we switch on <?ab? For a small inter-species repulsion, Qab = 
10~ 3 , all orbitals are still spread over two lattice sites and arrange in a staggered conformation 
as shown in Fig. 1C where there is half an A and half a B boson per site. For Qab = 0.1, 
all orbitals shrink to the size of a single site and the system is now in a Mott-insulator state 
where the A and B atoms arrange in alternating sites, see Fig. ID. In other words, phase 
separation on a length scale of one lattice size has happened. It is relevant to mention that an 
alternating pattern has been obtained by describing the bosonic mixture in the OL within an 
effective-spin Hamiltonian [9]. Finally, taking a much larger inter-species repulsion, Qab = 
10, a complete phase separation of the Mott-insulator state of the mixture emerges, see 
Fig. IE. It is instructive to examine the interface region of the two separated species, which 
shows no visible lowered density, contrary to the case of two phase-separated superfluids 
shown in Fig. IB. The reason is that the two-species Mott-insulator is incompressible. 

We have also investigated a mixture with the filling factors Ja,b = 1 and qa — Qb — 0.1. 
By minimizing Eab for Qab = we find that each boson resides in a different orbital located 
on a single lattice site. For qab = 10~ 3 we still have two Mott-insulator phases which sit 
atop one another, like in Fig. 3 A. For Qab = 0.1 we obtain that Eab is minimized by having 
two bosons per orbital. Hence, half the number of orbitals is required and, side by side, 
phase separation within neighboring sites occurs similarly to Fig. ID. Increasing the inter- 
species repulsion enforces particles of the same kind to share space and thereby enlarge the 
separation length-scale size. Finally, for a much larger inter-species repulsion, qab = 10, full 
phase separation occurs like in Fig. IE, but now with two bosons per orbital. 

Till now, we have described phase separations scenarios on length scales down to a single 
lattice site. As it turns, we cannot have phase separations on a shorter length scale if we 
are to stay in the 'lowest band' of the lattice. Let us now increase the interaction strengths 
such that 'higher bands' must be taken into account, see, e.g., [14], to properly describe the 
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ground state of the bosonic mixture. We consider the filling factors Ja,b = 1 and the intra- 
species parameters gA = 9b = 10. Obviously, for g AB — we have two strongly-interacting 
Mott-insulators which sit one atop the other. This is also the situation for g AB = 1, see 
Fig. 3A, and there is nothing in the appearance of these orbitals to tell us what is about to 
happen when g AB is increased. For g AB = 10 we do obtain phase separation on the shortest 
length scale, namely inside a single lattice site, see Fig. 3B. In each site, the two orbitals 
are pressed one against the other and thus narrow, and localize at the borders of the lattice 
site. In other words, the inter-species repulsion is sufficient to make the A and B bosons 
to separate in each site and form hetero-atom 'pairs'. What happens when g AB is further 
increased? For g AB = 30 we get an alternating arrangement of now homo- atom 'pairs', see 
Fig. 3C. Finally, for g AB = 40 the two strongly-interacting Mott-insulators fully separate in 
space as shown in Fig. 3D, which concludes our studies. 

Let us briefly summarize. We have shown that bosonic mixtures in optical lattices can 
undergo a wealth of phase separations on different length scales. The stronger the intra- 
species interactions are the shorter is this length scale. Specifically, (i) superfluids can only 
fully separate, i.e., occupy different regions of space; (ii) strongly-interacting superfluids 
can form alternating 'domains' whose sizes are a few lattice sites; (iii) Mott-insulators can 
undergo demixing where the atoms of each species reside in neighboring lattice sites; and, 
finally, (iv) strongly-interacting Mott-insulators can undergo demixing on the shortest length 
scale, namely, inside each lattice site. 



8 



REFERENCES 

[1] M. Greiner et al, Nature (London) 415, 39 (2002). 

[2] D. Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998). 

[3] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 (2006). 

[4] M. Boninsegni, Phys. Rev. Lett. 87, 087201 (2001). 

[5] E. Demler and F. Zhou, Phys. Rev. Lett. 88, 163001 (2002). 

[6] A. B. Kuklov and B. V. Svistunov, Phys. Rev. Lett. 90, 100401 (2003); A. Kuklov, N. 
Prokof'ev, and B. Svistunov, Phys. Rev. Lett. 92, 050402 (2004). 

[7] B. Paredes and J. I. Cirac, Phys. Rev. Lett. 90, 150402 (2003). 

[8] M. A. Cazalilla and A. F. Ho, Phys. Rev. Lett. 91, 150403 (2003). 

[9] E. Altman et al, New J. Phys. 5, 113 (2003). 

[10] K. Ziegler, Phys. Rev. A 68, 053602 (2003). 

[11] G.-P. Zheng, J.-Q. Liang, and W. M. Liu, Phys. Rev. A 71, 053608 (2005). 
[12] A. Isacsson et al, Phys. Rev. B 72, 184507 (2005). 
[13] M. P. A. Fisher et al, Phys. Rev. B 40, 546 (1989). 

[14] O. E. Alon, A. I. Streltsov, and L. S. Cederbaum Phys. Rev. Lett. 95, 030405 (2005). 

[15] T.-L. Ho and V. B. Shenoy, Phys. Rev. Lett. 77, 3276 (1996); B. D. Esry et al, ibid. 
78, 3594 (1997); H. Pu and N. P. Bigelow, ibid. 80, 1130 (1998). 

[16] L. S. Cederbaum and A. I. Streltsov, Phys. Lett. A 318, 564 (2003). 

[17] A. I. Streltsov, L. S. Cederbaum, and N. Moiseyev, Phys. Rev. A 70, 053607 (2004); O. 
E. Alon and L. S. Cederbaum, Phys. Rev. Lett. 95, 140402 (2005). 

[18] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

9 



FIGURES 



0.6 F 





A A A A 







/v A A A 





A A A A A A 




A A A A A A 



-2 















FIG. 1. (Color online) Phase separation scenarios of two superfluids for different intra- and 
inter-species interaction strengths. Shown are the orbitals </> and V' of the A and B species, re- 
spectively, and densities p (lower curves). Orbitals of A bosons are in black, gray, and blue, and 
of B bosons - in orange, magenta and red. Densities are in black (A bosons) and red (B bosons). 
The optical lattice is illustrated for guidance by the background sinusoidal curve. The index "i" 
enumerates lattice maxima. For convenience, the y/ 'rii / 'N ~A<t>i{x) an d Pa{x) = ^- Y^i° rb ^il'M^OI 2 
(scaled by |) are plotted, and similarly for the B bosons. All figures are for the filling factors 
Ia = Ib = \ in an optical lattice of Er = 20 recoil energies in depth and 16 lattice sites in length. 
The intra- and inter-species coupling parameters are: (A) gA = 9b = 10~ 3 , gAB = 10~ 4 . Orbitals 
of both species sit one on top of the other and cannot be distinguished; (B) gA = 9b = 10~ 3 , 
gAB = 1; (C) gA = gs = 0.1, gAB = 10~ 3 . Orbitals of both species sit one on top of the other with 
an 'offset' of one lattice site. For clarity of presentation B species orbitals are shifted upward and 
all orbitals are scaled by I; (D) g A = g B = 0.1, gAB = 0.1; (E) g A = g B = 0.1, gAB = 10. 
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FIG. 2. (Color online) Phase separations of two strongly-interacting superfluids. Same as in 
Fig. 1 except for the parameters: gA = 9b = 0.1 (all figures), and (A) qab = 0.1, Ja = Jb = |; (B) 
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10, f A = f B = h (C) <mb = 0.1, f A = f B = l (18 lattice sites); (D) 5ylfl = 10, = / B 



(18 lattice sites). 
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FIG. 3. (Color online) Phase separations of two strongly-interacting Mott-insulators. Same as 

in Fig. 1 except for the parameters: gA = 9b = 10, Ja = fs = 1 and no scaling of p (all figures), 
and (A) g AB = 1; (B) g AB = 10; (C) g AB = 30; (D) g AB = 40. 
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